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We have performed an angle resolved photoemission spectroscopy (ARPES) study of lightly-doped 
to underdoped YBa2CusO !/ (YBCO) untwinned single crystals and a core-level x-ray photoemission 
spectroscopy (XPS) study of YBCO single and polycrystals. In the zone diagonal (nodal) direction, 
dispersive quasi-particle (QP) features crossing the Fermi level were observed down to the hole 
concentration of ~ 4%, which explains the metallic transport of the lightly-doped YBCO. The 
chemical potential shift estimated from XPS was more rapid than in the Bi2212 cuprates. Upon 
hole doping, very rapid spectral weight transfer from high binding energies to the QP feature, even 
faster than L^-xSr^CuOi, was observed. 

PACS numbers: 74.25. Jb, 74.72.Bk, 79.60.-i 



INTRODUCTION 

The doping evolution of the electronic properties of 
lightly-doped cuprates have attracted considerable inter- 
est because of their fundamental importance in the con- 
text of both physics near Mott transition and mecha- 
nism of superconductivity.. In the cuprates, in particu- 
lar, the phase competition and the nature of the com- 
peting phases in the underdoped region have provided 
major experimental and theoretical challenges. Angle- 
resolved photoemission spectroscopy (ARPES) studies 
of lightly-doped La2- a; Sr :E Cu04 (LSCO) [l| have shown 
that a sharp dispersive quasi-particle (QP) peak crosses 
the Fermi level (Ep) in the nodal k = (0,0)-(7r,7r) direc- 
tion, which explains their metallic behavior |2[ and that 
the other part of the Fermi surface around (n,0) is gapped 
or pseudo-gapped, resulting in an "arc" of the Fermi sur- 
face around the nodal direction. A remnant of the lower 
Hubbard band (LHB) (or the polaronic side band 
persists at ~ —(0.4-0.6) eV as the QP peak appears near 
Ep. As holes are further doped, spectral weight is trans- 
ferred from the "LHB" to the QP peak 0). ARPES stud- 
ies on lightly-doped Caa^Na^CuOoCla (Na-CCOC) @ 
and Bi2Sr2- a; La a ;Cu06+5 (Bi2201) [a], on the other hand, 
have shown quite different behaviors: With hole dop- 
ing, the "LHB" feature is shifted toward Ep and merges 
into the QP band crossing Ep, making the LHB-to-QP 
spectral weight transfer less obvious. These observations 
are consistent with the pinning of the chemical potential 
in LSCO and the rapid chemical potential shift in Na- 
CCOC, Bi2201 and Bi2212 as implied by core-level x-ray 
photoemission spectroscopy (XPS) 0j 043] - 

Recent in-plane transport studies of YBa2Cu30j, 



(YBCO) have revealed their metallic behavior (dp/dT > 
0) at high temperatures even in the lightly-doped "in- 
sulating" region as in the case of LSCO M. The E F 
crossing of QP peak in the nodal direction has been ob- 
served in lightly-doped [10j as well as in underdoped to 
overdoped YBCO It is therefore important to see 

whether the electronic structure of lightly-doped YBCO 
is similar to that of LSCO or to that of Na-CCOC, 
Bi2201, and Bi2212. ARPES studies of YBCO have 
not been advanced compared to those of other high-T c 
cuprates because intense surface-state signals near Ep 
mask the bulk electronic states around the X [=(-7T,0)] 
and Y [=(0,7r)] points in the Brillouin zone. Lu et 



al. [12| have made high-resolution ARPES measurments 



on YBa2Cu3 06.993 and have succeeded in distinguishing 
bulk electronic states from the surface states around X 
and Y by aging the sample surfaces and deduced the 
superconducting gap. They have reported that the su- 
perconducting peak is accompanied by a peak-dip hump 
structure, similar to Bi2212. Yoshida et al. [10( and 



Borisenko et al. ll| have reported a bilayer splitting in 
lightly-doped and optimally-doped YBCO, respectively. 
In a recent ARPES study of optimally-doped YBCO by 
Nakayama et al. [3], bulk and surface Fermi surfaces 
have been separated. More recently, Hossain et al. re- 
ported the doping evolution of the electronic structure 
of the topmost Cu02 plane from overdoped to the un- 
derdoped region 14j. The recent Hall resistance mea- 
surment of underdoped YBCO revealed the existence of 
an small electron pocket [la,l3> although in the ARPES 
studies of other high-T c cuprates, only a Fermi arc around 
the nodal direction has been observed. Therfore, it is 
important to clarify whether an small electron pocket is 
observed in ARPES studies of underdoped YBCO or not. 
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TABLE I: Oxygen content y, hole concentration 8 in the 
CUO2 plane estimated from the electrical resistivity and ther- 
mopower and T c of the YBCO samples studied in the 
present work. 



y 





1c (K) 


6.20 


0.01 




6.25 


0.02 




6.28 


0.03 




6.30 


0.04 




6.35 


0.05 




6.40 


0.07 




6.45 


0.08 


20 


6.50 


0.09 


35 


6.60 


0.12 


57 



In this paper, ARPES and core-level XPS studies of 
lightly-doped to underdoped YBCO are presented to 
clarify the above issues and to gain further insight into 
the electronic structure of YBCO. 



EXPERIMENT 

Untwinned single crystals of YBa 2 Cu30 y (y=6.25, 
6.28, 6.30, 6.35, 6.40, 6.45, 6.60) were grown by the flux 
method as described in Ref. [13]. In the XPS study, 
we measured polycrystals of YBa2Cu30 y (y =6.20, 6.40, 
6.50, 6.60), too. In YBCO, how the in-plane hole concen- 
tration 6 changes with oxygen content y has been con- 
troversial. We have adopted the estimates from electri- 
cal resistivity and thermopower (Tsj . y, S and T c of the 
measured samples are listed in Table |U The ARPES 
measurements were carried out at BL10.0.1 of Advanced 
Light Source (ALS) using incident photons of 55 eV and 
a SCIENTA R4000 analyzer, at BL28 of Photon Fac- 
tory (PF) using incident photons of 65 eV, a SCIENTA 
SES2002 analyzer, and at BL5-4 of Stanford Synchrotron 
Radiation Laboratory (SSRL) using photons of 28 eV 
and a SCIENTA SES-200 analyzer. An R-Dec Co. Ltd. 
i GONIO LT goniometer ^ was used at PF. The total 
energy and momentum resolution was about 20 meV and 
0.027T in units of 1/a (a = 3.86 A is the lattice constant), 
respectively. The samples were cleaved in situ in an ultra 
high vaccum of 10" 11 Torr and cooled to ~ 10 K during 
the measurement. The XPS measurements were carried 
out using the Mg Ka line (hv = 1253.6 eV) and a VSW 
EAC 125 hemispherical analyzer for sintered polycrys- 
talline samples. Single cyrstals were measured using the 
Mg Ka line (hv = 1253.6 eV) and a SCIENTA SES-100 
analyzer. The energy resolution was about 0.8 eV, which 
was largely due to the width of the photon source. Nev- 
ertheless we could determine the core-level shifts with an 
accuracy of about ±50 meV because most of the spectral 
line shapes were identical between different compositions. 
Details of the core-level XPS studies are described else- 



where '8, 201. 



RESULTS AND DISCUSSION 

Figure [1] shows energy distribution curves (EDC's) 
along high symmetry lines in the second Brillouin zone 
for the y=6.35 and 6.60 samples. Two dispersive features 
are observed and both of them become more pronounced 
in going from y=6.35 to 6.60, reflecting the increased 
hole doping. One (open circles) is around —0.2 eV at T, 
dispersing upward and crossing Ep around 60 % of the 
rX and TY lines. In optimally-doped YBCO, this feature 
was assigned to surface states in R ef. J 1 21] and assigned to 
surface antibonding band in Ref. [13| The other feature 
(open triangles) was considered to be the surface bonding 



band according to Ref. [13| ■ The dispersion of the surface 
bonding band is unclear around T because the broad and 
strong antibonding band masks this feature. It should 
be noted that the energy position of the surface bonding 
band showed clear anisotropy, —0.18 eV at X and —0.14 
eV at Y. The ID like chain-derived states which was ob- 
served in the spectra of optimally doped YBCO [IH is 
not observed here, perhaps due to the oxygen deficiency 
in the CuO chain and also due to the difference in the 
photon energy and the polarization of the incident light, 
i.e., difference in the matrix elements. 

Figure [UJa) shows the plot of spectral weight in the 
second Brillouin zone integrated within 40 meV of Ep, 
which reflects the Fermi surface, for y=6.45. Not only a 
large hole-like Fermi surface centered at the S point, but 
also an apparently electron-like Fermi surface centered at 
T are observed, indicating the bilayer splitting loL 11|. 
Open circles in Fig.[2ta) are the hp positions of the Fermi 
surfaces determined by the peak positions of momentum 
distribution curves (MDC's) at Ep. In Fig.^b), the k F 
positions thus determined are plotted for various compo- 
sitions in the same panel. The shape of the Fermi surface 
does not vary appreciably in this doping range, suggest- 
ing that these Fermi surfaces are not derived from bulk. 

Because the cleaved surface of YBCO tend to be over- 
doped 11, llij . we compare the experiment with the 



Fermi surfaces predicted by local-density-approximation 
(LDA) calculation for y=7 [5l| in Fig.^b). In the LDA 
calculation, the Fermi surface has four sheets derived 
from the Cu02 planes and the CuO chain. However, be- 
cause we did not observe the CuO chain-derived bands, 
we have plotted only the Cu02 plane-derived sheets, that 
is, those derived from the bonding and antibonding bands 
of the Cu0 2 bilayer. The Fermi surfaces for two values 
of fc z =0 and ir/c are indicated because for the optimally- 
doped YBCO the CuO chain mediates orbital overlap 
along the c-axis and therefore strong k z dispersions are 
predicted. The experimental results are more similar to 
the fc z =0 cross-section than to the k z —n/c one. The 
shape of the Fermi surface of the bonding band is more 
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similar to that of Bi2212 than to that of LSCO. According 
to the tight binding model including hopping parameters 
t, t', and t", this suggests a large next-nearest-neighbour 
hopping \t'\ in YBCO @. 

In Fig. |2Ja), we summarize the surface band disper- 
sions of YBCO for various y's obtained from EDC's. 
Here, the surface bonding band for the ?/=6.28 sample 
was taken in a different experimental geometry {E _L b) 
and photon energy {hv — 65 eV) than the rest of the data 
(E || XY), and was clearly resolved even around the T 
point. Although the spectral line shape and intensities 
dramatically changed with hole doping (see Fig. [T]), the 
band dispersions did not change appreciably in this dop- 
ing range. The lower bound for the bilayer splitting at 
X and Y is estimated to be 180 meV and 140 meV, re- 
spectively, because the surface antibonding band should 
be above E-p at these momenta. Figure [3jb) shows the 



LDA band dispersions of YBCO (y=7) for k z =0 [21 1. 
Red curves are the bonding and antibonding bands of the 
Cu02 bilayer. The experimental surface band dispersions 
are qualitatively reproduced by the LDA band structure 
except for the E-p crossing points along the symmetry 
lines and the overall band narrowing (by a factor of ~ 0.5) 
in experiment. This similarity indicates that the origin 
of the surface state is the overdoped cleaved surface. [l3j 
In the LDA calculation, the bonding- antibonding split- 
ting of the Cu02 bilayers at X and Y was 610 meV 
and 560 meV, respectively, which is much larger than 
the experimental value of > 180 meV and > 140 meV. 
The reason for this dicrepancy is that LDA calculation 
neglects electron correlation and generally overestimates 
band dispersions. Indeed, overdoped Bi2212 has shown 
a bilayer splitting of 88 meV [2J|, much smaller than 
300 meV predicted by LDA calculation [2l|,[25j]. Bilayer 
Hubbard model calculations for two coupled Cu02 layers, 
which explicitly include the on-site Coulomb repulsion, 
predicted a maximum splitting of 40 meV [2(| . 



=6.35 5-0.05|y=6.60 8-0.12] 



ri , Ba 2 Cu 3 O s ,hv=55 eV | 





-0.4 -0.8 -0.4 0-0 
Energy relative to E F (eV) 



FIG. 1: (Color online) EDC's of lightly-doped to underdoped 
YBCO. Two features in the spectra are surface bonding state 
(open triangles) and surface antibonding state (open circles). 
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FIG. 2: (Color online) (a) Spectral weight mapping at E-p for 
a 40 meV integration window. Open circles are the A?f posi- 
tions determined from MDC's and the white arrow denotes 
the polarization of the incident light, (b) Experimentally de- 
termined &f positions in each doping level and LDA Fermi 
surface for k z —0 and fc z =7r/c related to the Cu02 plane. 





FIG. 3: (Color online) Band dispersions in YBCO obtained 
from EDC's (a) compared with the LDA band dispersions for 
YBCO (y=7) at fc z =0 [HI (b). In (b), red curves are the 
antibonding and bonding bands of the Cu02 planes. 



Figure 4(a)-(c) shows spectral weight mapping at E-p 
with a 40 meV integration window in the first Brillouin 
zone for the y=6.35, 6.45 and 6.60 samples. A Fermi 
"arc" around the nodal direction is observed as in the 
LSCO case, which is completely different from Fig. 2(a). 
Figure 2(a) does not show a Fermi arc but a hole-like 
and an electron-like Fermi surface arising from the sur- 
face state. Figure 4(e) shows the intensity plot in the 
E-k plane in the nodal direction for the y=6.35 sam- 
ple together with the energy band dispersion determined 
from MDC peaks. The dispersion of this band shows 
a kink structure around ~ 50 meV like in many other 
high-T c cuprates, while the surface bonding band is re- 
ported to show a straight dispersion [13| . In addition, 
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the &f positions in the nodal direction in the first Bril- 
louin zone are different from those in the second Brillouin 
zone. These findings suggest that the Fermi arc observed 
in the first Brillouin zone is distinct from the Fermi sur- 
faces observed in the second BZ and we attribute it to 
a bulk-derived feature. To investigate the doping evolu- 
tion of the Fermi arc, the kp positions dctcrmicnd by the 
peak positions of the MDC's are plotted by filled circles 
in Fig. 4(a)-(c). The fep positions for y—6.35, 6.45, and 
6.60 are overlaid in Fig. 4(d). As holes are doped, the 
hole-like Fermi surface centered at S slightly expands in 
the nodal region, which is also consistent with the as- 
sumption that these Fermi surfaces are bulk derived. 
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FIG. 4: (Color online) (a)-(c) Spectral weight mapping at Ef 
for a 40 meV integration window in the first Brillouin zone. 
Filled circles are the fcp positions determined from MDC's 
and the white arrow denotes the polarization of the incident 
light, (d) fcp positions for i/=6.35, 6.45, and 6.60 determined 
from MDC's. (e) Intensity plot in the E-k plane in the nodal 
direction for the y=6.35 sample. Filled circles denote the 
energy band dispersion determined by fitting the MDC's to 
Lorentzians. 



Next, we present the result of the measurements of 
XPS core levels in YBCO, and deduce the chemical po- 
tential shift as a function of doping. As shown in the 
inset of Fig. El the Ols and Y 3c? core levels show the 
same shifts. As for the Ba Ad core level, the peak be- 
comes broader with hole doping, perhaps due to some 



surface effects, but apart from this the Ba 4c? peak po- 
sition follow the shifts of the Ols and Y 3 c? core levels, 
too. The line shape of the Cu 2p peak was not identical 
between the different compositions but its peak position 
moved in the opposite direction to those of the O Is, 
Y 3d and Ba Ad core levels. The different behavior of 
the Cu 2p core level can be attributed to the creation 
and the growth of "Cu 3+ " signals with hole doping on 
the higher binding energy side of the Cu 2+ main compo- 
nent [4|. The inset of Fig. [5] shows the binding energy 
shifts referenced to y=6.2 of each core level thus esti- 
mated. As in the previous studies [2(j, we consider that 
the parallel shifts of the O Is, Y 3d and Ba 4c! core levels 
reflect the chemical potential shift A/i. Recently, Maiti 
et al. reported the shifts of the core levels in YBCO us- 
ing hard x-ray photoemission spectroscopy [271 ] ■ Our Y 
3c? and O Is shifts are consistent with their results ex- 
cept for their lowest doped sample (y — 6.15). In the 
main panel of Fig. [5j therefore, A/i has been obtained by 
taking the average of the Ols and Y 3d core- level shifts 
as shown in the figure. In the same panel, the chemical 
potential shift in other high T c cuprates are plotted for 
comparison 0, 0, 1, Hi- Theoretically, if the Ion g-range 
antiferromagnetic order of the parent Mott insulator dis- 
appears for a very small hole concentration, the charge 
suceptibility dn/dfi diverges in the limit 5 — > [29j and 
hence the chemical potential is suppressed. In LSCO 
and Bi2212, the antiferromagnetic order quickly disap- 
pears for a small amount of hole doping and indeed the 
shift is suppressed for small S. In YBCO and NCCO, on 
the other hand, the antiferromagnetic order perisits up 
to (5 ~ 0.07 and <5 ~ 0.14, respectivly, and thefore it is 
reasonable that little or no suppression of the shift was 
observed for small <5. 

According to t-t'-t"-J model calculations [30], the 
larger the \t'\ becomes, the stronger the chemical po- 
tential shift becomes. The observed differences between 
Bi2212 and LSCO, such as the different dispersion width 
of the LHB from k ~ (tt/2, tt/2) to ~ (tt, 0), the different 
shape of the Fermi surface as well as the different behav- 
ior of the chemical potential shift, can be explained by a 
larger \t'\ value in Bi2212 than that in LSCO [3(J ffl^. 
A/z in YBCO is the largest among those compounds, 
which means that the |t'| of YBCO may be the largest, 
consistent with the conclusion from the ARPES results 
described above. The suppression of the shift in the un- 
derdoped LSCO is too strong compared to the t-t'-t"-J 
model calculation and has been attributed to a fluctua- 
tion of charge stripes [3]. The dynamical stripes seen in 
LSCO change their periodicity linearly with doping (32^ . 
As a result, the local charge density is nearly unchanged 
so that the hole rich region expands and the hole poor 
region shrinks as the average hole concentration is in- 
creased. Under such a circumstance, the chemical poten- 
tial would not move with doping. If the periodicity of 
the stripes remains unchanged with doping, on the other 
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hand, the chemical potential is expected to move in a 
rigid-band-like manner. In this regard, the observations 
of in-plane anisotropic magnetoresistance [33| and elec- 
trical resistivity Q in untwinned single crystals of YBCO 
are most consistent with the idea that charge stripes also 
exist in lightly-doped YBCO. The presence of the stripes 
and the strong chemical potential shift can be reconciled 
if the stripes in the lightly-doped YBCO are of the type 
of fixed periodicity. 
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FIG. 5: (Color online) Chemical potential shift in YBCO com- 
pared with those in other high T c superconductors pi. 17l.lsl.l2sl]. 
The antiferromagnetic region is displayed by shaded region for 
YBCO and NCCO. Inset shows the shift of each core level for 
YBCO as a function of hole concentration S. 



To investigate the doping dependence of the shape of 
the Fermi surface and the QP band quantitatively, we 
have fitted the fcp positions shown in Fig. 4(d) to the 
two-dimensional t-t'-t" tight binding model 

£fc — M = — 2t[cos{k x a) + cos^a)] — 4i'cos(/c :! ;a)cos(£;j,a) 
— 2t" [cos(2k x a) + cos(2k y a)] + e n 

Here, en is the position of the band center relative to the 
chemical potential. We have assumed that t"/t' = — 1/2 
as before @, Ho|, 34 1, and regarded —t'/t and — eo/t 
as fitting parameters. Figure 6(a) shows the shape of 
the Fermi surface from the tight-binding fit. One can 
see the squarelike Fermi surface which suggests large 
—t'/t pH - Figure 6(b) shows fitted tight-binding param- 
eters as function of hole concentration. As in the case 
of LSCO and Bi2201, —t'/t shows weak doping depen- 
dence @, [34|. However, the —t'/t value is much larger in 
YBCO (~ 0.4) than that in LSCO (~ 0.2) and Bi2201 (~ 
0.2), consistent with the XPS results. If the chemical po- 
tential moves in a rigid-band-like manner, A/i/t = eo/t, 
as is not the case in YBCO and LSCO. The doping depen- 
dence of eo/t in YBCO is similar to that in LSCO while 
—Afi/t increases faster than that in LSCO. This differ- 
ence may be attributed to the different type of charge 
stripes as described above. 
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FIG. 6: (Color online) Doping dependence of the electronic 
structure in YBCO. (a) Doping dependence of the Fermi sur- 
face shape, (b) Doping dependence of the fitted tight-binding 
parameters. 
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FIG. 7: (Color online) ARPES spectra taken at 10 K for 
lightly-doped YBa2Cu30 H along the nodal direction in the 
first Brillouin zone. S indicates hole concentration per CUO2 
plane, (a) Top panel shows EDC's and bottom panel shows 
intensity plot in the E-k plane, (b) EDC's at kp for y=6.25, 
6.30, 6.35, and 6.40 samples, (c) fep positions for various 
doping levels in YBCO, LSCO Qj, Na-CCOC 0, B12201 @|, 
and Bi2212 H. 



Figure EJa) shows EDC's (top) and corresponding in- 
tensity plots in E-k space (middle) in the nodal (0,0)- 
(tt,tt) direction of the y=6.30 sample. As in the case of 
LSCO, a dispersive QP feature crossing E-p is seen for 
the lightly-doped y >6.30 (<5=0.04) sample and its inten- 
sity increased with doping as shown in Fig. [Tfb), giving 
rise to the metallic transport at high temperatures [§]. 
Figure [7£c) shows the k-p positions along the nodal di- 
rection for various doping levels in YBCO, LSCO, Na- 
CCOC, Bi2201, and Bi2212. In the antiferromagnetic 
(AF) insulating state of the undoped compound, the 
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valence-band maximum along the nodal direction should 
be located at the AF zone boundary, that is, exactly at 
(7r/2,7r/2) due to the band folding in the AF state. In 
Na-CCOC and the Bi-based compounds, indeed, with de- 
creasing hole concentration the kp position moves toward 
(7r/2, 7r/2). In LSCO, on the other hand, kp extrapolates 
to ~(0.447r, 0.44tt) and not to (7r/2,7r/2) [34]. In LSCO, 
in addition to the QP feature crossing Ep, the broad 
"LHB" or the polaronic feature at ~ -(0.4-0.6) eV is rec- 
ognized for low hole concentrations x < 0.05. Such an 
electronic structure has been regarded as the presence of 
the coherent and incoherent parts of the spectral func- 
tion influenced by electron-electron and electron-phonon 
interacitons Q and/or due to a nano-scale electronic in- 
homogeneity such as charge stripes [H]. Since such an 
inhomogeneity is also expected to exist in lightly-doped 
YBCO ESI. the fact that fc F in YBCO did not extrapo- 
late to (7r/2, 7r/2) suggests that QP band is also separated 
from the "LHB" in YBCO. Unfortunately, such a rem- 
nant of the "LHB" is not seen in YBCO at least down 
to ~ —0.80 eV, partly because it is masked by the strong 
overcapping tail of the O 2p band and partly because 
spectral weight transfer from the LHB to the QP band 
occurs very quickly with hole doping. 

CONCLUSION 

We have observed the band dispersions and Fermi sur- 
faces of lightly-doped to undcrdoped YBCO by ARPES 
using untwinned single crystals. The shape of the Fermi 
surface and the large chemical potential shift indicated 
large like Bi2212, Bi2202, and Na-CCOC. On the 
other hand, we observed a rather clear dispersive QP 
feature crossing E-p along the nodal direction, revealing 
a similarity to LSCO. Although the coexistence of the 
remnant LHB could not be confirmed unlike the other 
cuprates, we conclude that very rapid spectral weight 
transfer ocurrs from the LHB to the QP band upon hole 
doping in the lightly-doped to underdoped YBCO. We 
have also studied the chemical potential shift with hole 
doping by the XPS measurements of single and poly- 
crystals. The shift was as rapid as those in Na-CCOC, 
Bi2201, and Bi2212 in spite of the fact that charge stripes 
are observed in lightly-doped YBCO like LSCO. We at- 
tribute the different chemical potential shifts between 
YBCO and LSCO to different types of charge stripes. 
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